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ABSTRACT 



Aims. We compare the predictions of evolutionary models for early-type stars with atmospheric parameters, projected 
rotational velocities and nitrogen abundances estimated for a sample of Be-type stars. Our targets are located in 4 
fields centred on the Large Magellanic Cloud cluster: NGC2004 and the Nil region as well as the Small Magellanic 
Cloud clusters: NGC 330 and NGC 346. 

Methods. Atmospheric parameters and photospheric abundances have been determined using the non-LTE atmosphere 
code tlusty. Effective temperature estimates were deduced using three different methodologies depending on the 
spectral features observed; in general they were found to yield consistent estimates. Gravities were deduced from 
Balmer line profiles and microturbulences from the Si ill spectrum. Additionally the contributions of continuum emission 
from circumstellar discs were estimated. Given its importance in constraining stellar evolutionary models, nitrogen 
abundances (or upper limits) were deduced for all the stars analysed. 

Results. Our nitrogen abundances are inconsistent with those predicted for targets spending most of their main sequence 
life rotating near to the critical velocity. This is consistent with the results we obtain from modelling the inferred 
rotational velocity distribution of our sample and of other investigators. We consider a number of possibilities to 
explain the nitrogen abundances and rotational velocities of our Be-type sample. 
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Magellanic Clouds 



1. Introduction 



stars: emission-line, Be - stars: atmospheres - stars: rotation - stars: evolution - galaxies: 



The classical Be-type star is a B-type star that shows or 
has previously shown prominent emission features in its 
Balmer line spectrum, indicating the presence of a geo- 
metr ic ally f lattened, circumstellar disc ()Quirrenbach et ail 
11994 119971 ). Although similar phenomena can be ob- 
served in some late O- and earl y A-type stars , it is more 
prevalent i n B spectral types (jConti fe Leed 119741 - O- 



type stars; iHerbig] 119601 - A a nd B-type emi s sion stars). 
IJaschek fc Jaschekl (jl983h and IZorec fc Briotl (|l997l ) have 
estimated that the frequency of Galactic stars that show 
the Be phenomenon is approximately 17%, with the high- 
est fraction around spectral types Ble - B2e. Emission is 
also often seen in the helium and iron spectra, with sili- 
con and magnesium emis sion lines being seen in some stars 
([Porter fc Riv inius 2003ft , Be-type stars can show a variety 
of spectral characteristics and some can lose their emission 
line spectrum and transition to a normal B spectral type 
(|Claretil2003h . 

Although the mechanism of the Be phenomena is still 
unclear, such stars are beli eved to rotate with veloci- 
ties > 40% (ICranmerl l2005h up to - 90% of the crit- 
ical lim it, where the centrifu gal force balances that of 
gravity (|Townsend et al.ll2004l ). The Magellanic Clouds are 



ideal stellar laboratories to investigate the effects of rapid 
rotation. These lower metallicity regimes generally pro- 
duce stars with weaker stellar winds and this reduces the 
loss of angular momentum, leading to a larger fraction of 
rapidly rotating stars. Additionally at lower metallicity, 
stars of a given mass are more compact again resulting 
in larger rota t ional veloc ities. Therefore, as discu ssed by 
IMaeder et all (|l999tl and iMaeder fc Mevnetl (|2000D . lower 
metallicity re gimes should support a larger f raction of Bc- 
type objects. IWisniewski fc BiorkmanI (|2006l ) found an in- 
crease in the proportion of Be to B-type stars with de- 
creasing metallicity, whilst iMartavan et all (|2006ll2007l ) ob- 
served a large sample of B and Be-type populations in 
the Magellanic Clouds and compared them with a Galactic 
sample, concluding: "the lower the metallicity is, the higher 
the rotati onal velocities are" . Th is conclusion was sup- 
ported bv lHunter et al.l I 



* Based on observations at the European Southern 
Observatory in programmes 171.0237 and 073.0234 



but Keller et al.l ([1999D ob- 
served no difference in their fraction of Be-type stars for 
different metallicity regimes as part of a photometric based 
study. 

The ev olutionary s t atus o f Be-type stars has been dis- 
cussed bv IZorec et al.l (|2005l ). By determining the stellar 
mass and age of 97 Be-type objects they reached two im- 
portant conclusions: Be-type stars are not constrained to a 
single stage of the main-sequence but are present over the 
whole of this evolutionary phase, and more massive stars 
(M > 12 Mq) present the Be phenomena typically earlier in 
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their evolution than less massive stars. They conclude that 
this arises from the high mass-loss rates of more massive 
objects that can reduce their speed of rotation, and also 
the time scales involved for transferring angular momen- 
tum from the core to the surface, which is longer for lower 
mass stars. 

The sample of Be-type stars considered here have been 
briefl y discussed as part of a larger B-type star survey 
by iHunter et all ([20081. l2009al ). wh o undertook a qualita- 
tive c omparison with the results of Martavan et al. (2006, 
12007ft . In general good agreement was found e.g. more Be- 
type stars are present in the SMC than LMC and their 
projected rotational velocities are higher than B-type stars. 
This paper develops and extends the work of Hunter et 
al. by estimating the atmospheric parameters and chemical 
compositions for this sample. 

Nitrogen can be utilised to investigate the degree of mix- 
ing of nucleosynthetically processed material from the stel- 
lar core into the atmosphere. Such material shows large en- 
hancements in its nitrogen content, together with smaller 
changes in the amounts of carbon, oxygen and sodium, 
other lighter element s such as boron, lithium and beryllium 
can also be affected (iHeeer fe L angcr 200j| iPrzvbilla et al.l 
l2010t iBrott et al.ll2011a|) . Hence nitrogen is ideal for testing 
stellar evolutionary models that include the effects of rota- 
tional mixing and here we estimate nitrogen abundances 
within a sample of 61 Be-type stars from observation fields 
centred on four open clusters, two in the LMC and two in 
the SMC. In Sect. [21 we outline the observational data col- 
lected as part of th e VLT-FLAMES survey of massive stars 
(|Evans et al.1 12005^. Sect.[3]d iscusses the procedures used to 
obtain atmospheric parameters and abundance estimates, 
together with a description of the methodology used to cor- 
rect our data for contamination due to a circumstellar disc. 
Our results are presented and discussed in Sect. 0] and [5] 



2. Observations 

The spectroscopic data were obtained during a 
European Southern Observatory (ESO) large programme 
(|Evans et al.l I2005D using the Fibre Large Array Multi- 
Element Spectrograph (FLAMES) on the Very Large 
Telescope. Observations were obtained for seven clusters 
within three different mctallicity regimes (Galactic, LMC 
and SMC). A summary of the spectroscopic configurations 
is presented in Table [TJ with full details of target selection, 
data reduct i on and obs ervational details being given in 
lEvans et al.l (|2005l 120061 ). 



2.1. Selection criteria 

The focus o f pre vi ous p apers has been O-type 
(|Mokiem et al.l [20051 [20061 l2007h. narrow lined B- 



tvpe stars ([Trundle et alF 2007 HHunter et a 



I2007h or 

fast rotating B-type stars (see Hunter et al.l 2009aL and 
references therein). This paper will provide an analysis 
of the Be-type stars that were previously omitted from 
these investigations. For the Galactic clusters (NGC3293, 
NGC4755 and NGC6611), no such objects were observed 
and therefore all our targets lie within the Magellanic 
Clouds. The sample is summarised in Table [2j the spectral 
type and projected rotati onal velocity (v sin i ) being taken 
from [Hunter et al. (2008) unless otherwise stated. 



2.2. Equivalent width measurements 



Equivalent widths (EWs) for metal absorption lines in 
narrow lined spectra can normally be reliably measured 
by fitting Ga ussian profiles to th e observed features (see, 
for example FTrundle et all |2007|) . However as noted in 
IHunter et al.l (|2009al) . for stars with projected rotational 
velocities, usini> SOkms" 1 , the dominant broadening 
mechanism is rotation. Hence it was necessary to utilise ro- 
tationally broadened profiles in order to reliably estimate 
the equivalent widths of our sample of Be-type stars. 

The high rotational velocities of Be-type stars lead to 
additional c omplications when an alysi ng the strength of 
metal lines. IVrancken et al.l (|l997l ) and iPorter fc Riviniusl 
(2003) discuss these and in particular that the continuum 
can mask the outer wings of features present in the spec- 
tra of such objects. This can affect the estimation of the 
continuum level, line width and the line strength, normally 
resulting in a systematic u nderestimation of th ese quan- 
tities. To address this issue, IHunter et al.l ([20081 ) deduced 
projected rotational velocities from the helium spectrum 
which is less susceptible to such effects due to its intrinsic 
strength. As noted in Table [2[ six objects have had their 
rotational velocity estimates re-determined and these have 
been used throughout the current analysis. This was be- 
cause rotati onally broadened p rofiles generated using the 
estimates of IHunter et al. (2008) did not reliably reproduce 
the observed metal lines. Typically the difference was ^30% 
of the original value. 

Equivalent width measurements of metal lines were 
then undertaken u sing the methodology discussed by 
IHunter et al.l (|2007l ). Briefly, single or multiple rotational 
broadened profiles generated with the projected rotational 
velocities listed in Table [2] were fitted to the observed spec- 
tra using a least-squares technique. The significant broad- 
ening can cause blending with other lines, which can also 
affect the equivalent width measurements. The effects of 
such blending is illustrated for the star, NGC 346-008 (with 
vsvui^ 299 kms^ 1 ) in Fig. [T] for the Si in triplet near 
4560A. Lines were excluded from the subsequent analysis if 
rotationally broadened profiles that had been constrained 
to the appropriate wavelength centroids did not yield a sat- 
isfactory fit. 

Table [3] lists the equivalent width estimates for selected 
lines of nitrogen, magnesium and silicon. NGC 346-096 has 
been included in this table due to the measurement of N II 
line, but the lack of reliable Siin/Mgii spectra meant it 
was excluded from further analysis. One of the primary 
aims of this paper is to determine reliable nitrogen abun- 
dances for our sample. Hence for targets where no Nil line 
at 3995A could be distinguished, upper limits to this equiv- 
alent width were estimated. The methodolog y adopted has 
been m odified from the approach discussed in IHunter et al.l 
(|2007l ). as follows. On a star by star basis, a rotationally 
broadened profile was added to the observed spectrum and 
its strength varied until it was visible within the spectral 
noise. The line was then measured and the resulting equiv- 
alent width (rounded to the nearest 5 mA) was taken as a 
conservative upper limit. 
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Table 1. Wavelength coverage, mean FWHM of the arc lines and effec tive resolving power, R for each Giraffe central 
wavelength setting, A c towards NGC346, taken from lEvans et all (|2006l ). 



Setting 


Ac 


Wavelength Coverage 


FWHM 


R 




(A) 


(A) 


(A) 


(pixels) 




HR02 


3958 


3954 - 4051 


0.18 


3.7 


22 000 


HR03 


4124 


4032 - 4203 


0.15 


3.4 


27 500 


HR04 


4297 


4187 - 4394 


0.19 


3.5 


22 600 


HR05 


4471 


4340 - 4587 


0.23 


3.7 


19 450 


HR06 


4656 


4537 - 4760 


0.20 


3.5 


23 300 


HR14 


6515 


6308 - 6701 


0.39 


3.9 


16 700 



Table 2. List of Be-type stars observed in the Magellanic Cloud clusters NGC 346, NGC 330, NGC 2004 and the Nil 
region. 



Star 


Spectral 


S/N 


v sin i 


Note 


Star 


Spectral 


S/N 


v sini 


Note 




Type 




(kms -1 ) 






Type 




(kms -1 ) 




NGC 346-004 


Be(Bl:) 


307 


266 




NGC 330-065 


Bl-3e 


39 


284 


1 


NGC 346-008 


Ble 


187 


299 




NGC 330-068 


B1.5e 


30 


23 


1 


NGC 346-009* 


BOe 


114 


181 




NGC 330-069 


B3 Hie 


40 


193 




NGC 346-023* 


B0.2e 


130 


50 




NGC 330-070 


B0.5e 


39 


348 


1 


NGC 346-024 


B2e:shell 


160 


190 




NGC 330-076 


B3e 


35 


62 


1 


NGC 346-036 


B0.5 Ve 


142 


287 




NGC 330-083 


B3 Hie 


38 


140 




NGC 346-041 


B2e 


130 


144 


1 


NGC 330-085 


B3:e 


38 


191 


1 


NGC 346-045 


B0.5 Vne 


135 


181 




NGC 330-087 


Be-Fe 


24 


214 


1 


NGC 346-048 


Be (B3 shell) 


121 


158 


2 


NGC 330-091 


BOe 


32 


272 


1 


NGC 346-061 


Bl-2e 


118 


336 


1 


NGC 330-096 


Bl-3e 


26 


175 


1 


NGC 346-064 


Bl-2e 


117 


108 


1 


NGC 330-100 


Be (BO-3) 


25 


373 


1 


NGC 346-065 


B3e 


105 


222 


1 


NGC 330-112 


Bl-3e 


24 


262 


1 


NGC 346-067 


Bl-2e 


105 


351 


1 












NGC 346-068 


B0 Ve 


109 


378 


1 


NGC 2004-023 


B2e 


169 


102 




NGC 346-069 


Bl-2e 


102 


186 




NGC 2004-025 


B2e 


153 


83 




NGC 346-072 


Bl-2e 


112 


102 




NGC 2004-027 


BOe 


148 


182 




NGC 346-073 


Bl-2e 


110 


190 




NGC 2004-035 


Ble 


125 


244 




NGC 346-076 


B2e 


106 


237 


1 


NGC 2004-039 


B1.5e 


106 


212 




NGC 346-089 


Bl-2e 


91 


79 




NGC 2004-048* 


B2.5e 


86 


244 




NGC 346-091 


Ble 


77 


49 




NGC 2004-056 


B1.5e 


74 


229 




NGC 346-095 


Bl-2e 


97 


227 


1 


NGC 2004-058 


09.5 Ve 


125 


195 




NGC 346-096 


Bl-2e 


100 


343 


1 


NGC 2004-067 


B1.5e 


56 


237 




NGC 346-110 


Bl-2e 


93 


243 


1 


NGC 2004-083 


B1.5:e 


73 


194 


1 












NGC 2004-089 


B2.5e 


67 


288 




NGC 330-025* 


B1.5e 


93 


250 




NGC 2004-092 


B2e 


61 


171 




NGC 330-029 


B0.2 Ve 


85 


209 




NGC 2004-096 


B1.5e 


63 


245 




NGC 330-031 


B0.5 Ve 


98 


178 




NGC 2004-115* 


B2e 


111 


15 


3 


NGC 330-034 


Bl-2e 


54 


231 


1 












NGC 330-044 


Bl-2e 


25 


184 


1 


N 11-073 


B0.5e 


143 


222 




NGC 330-050 


B3e 


67 


214 


1 


N 11-074 


B0.5e 


149 


102 




NGC 330-054 


B2e 


44 


147 


1 


N 11-078 


B2e 


153 


93 


2 


NGC 330-060 


B2.5e 


37 


88 


1 


N 11-081 


B0:ne 


135 


363 


1 


NGC 330-062 


B3e 


43 


241 


1 


N 11-090* 


B2e 


119 


118 




NGC 330-064 


B3:e 


23 


269 


1 













Notes. Spectral types and projected rotational velocities (wsini) are taken from iHunter et aD (|2008f ). unless stated otherwise. 

Value of v sin i has been redetermined (see Sect. 12.2]) ™ Omitted: No reliable measurements of the Mgll/Silll or Hell spectrum 
' 2 ' Omitted: The presence of a strong shell spectra Omitted: Not a classical Be-type star. Short period binary. 



SPEC (|Hubenvl [19881: iHubenv fc Landfl99l iHubenv et all 
I1998D . The grid used covers the effective temperature range 
12 000 - 35 000 K in steps of 1 500 K at its lowest resolution 
and 500 K at its highest resolution; surface gravity covers 
the range 4.5 dex down to the Eddington limit in steps of 
0.25 dex. The microturbulent velocities have discrete values 



3. Analysis 

3.1. Non-LTE atmosphere calculations 

The analysis in this paper follows that of previous work on 
the FLAMES dataset, making use of a B-type stellar non- 
LTE model atmosphere grid, computed using the tlusty 
model atmosphere codes and the line formation code SYN- 
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Table 3. Equivalent width measurements and associated uncertainty estimates for 32 Be-type stars. 



Star 


Nn 3995 


Sim 4552 


Sim 4567 


Sim 4574 


Si iv 4116 


Sin 4130 


Mgn 448: 




(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


NGC 346-004 


133±7 


211±10 


154±8 


84±5 


83±5 




61±4 


NGC 346-008 


65 ±7 


274±15 


216±12 


106±6 






112±6 


NGC 346-009 


42 ±4 


101±8 


81±6 


49±4 


163±13 




31±2 


NGC 346-023 


31±5 


74±13 


52±9 


34±6 


49±9 


24±4 




NGC 346-024 


71±14 


68±14 




50±10 




87±18 


135±28 


NCC! 346-036 

IN V_TV_,' J1U UOU 


68±14 


88±18 


117 ±22 


106±21 








NGC 346-045 


60±8 


93±11 


77±9 


52±6 




36±10 


52±6 


NGC 346-069 


<40 


99±18 


68±13 


41±14 




63±12 


52±14 


NGC 346-072 


<35 


48±10 


27±8 


16±9 






51±10 


NGC 346-073 


<40 


71±13 


63±14 


36±11 






57±8 


IN VjrVy OIU UOy 


<25 


70+7 












NGC 346-091 


<25 


75±11 


58±8 


46±5 






45±6 


NGC 346-096 


120±17 














NGC 330-025 


56±10 


171±10 


105±11 


70±10 






101±13 


NGC 330-029 


<50 


53±14 


30±11 


27±10 


38±8 






NGC 330-031 


<40 


145±30 


87±18 


47±10 


51±10 






NGC 330-069 


<90 


73±13 


46±13 








160±30 


NGC 2004-023 


39±8 


78±15 


74±15 


25±6 






61±12 


NGC 2004-025 


24±5 


74±15 


63±13 


32±6 






99±19 


NGC 2004-027 


54±3 


102±6 


76±5 


59±5 








NCC 9004-035 


75±15 


29±6 


57±11 


30±6 






90±17 


NGC 2004-039 


52±6 


133±16 


110±13 


49±7 






195±23 


NGC 2004-048 


57±8 


116±13 


108±12 


70±10 






202±22 


NGC 2004-056 


53±14 


129±30 










119±17 


NGC 2004-067 


83±15 


163±29 


135±24 


94±17 






135±24 


NGC 2004-089 


36±10 


110±14 


94±12 


60±11 






182±24 


NGC 2004-092 


<40 


104±14 


96±12 


69±22 






114±30 


NGC 2004-096 


54±16 


144±43 


101±30 


77±23 






82±25 


N 11-073 


32±9 


54±15 










52±14 


N 11-074 


30±6 


84±16 


92±18 


81±16 


58±12 






N 11-090 


19±5 


67±14 


68±14 


34±7 






52±11 



Notes. The Sim and Mgn absorption lines listed were necessary in the methodology used to obtain atmospheric parameters 
(Sec. GO]). Nil at 3995 A was used to compare our sample with evolutionary models (Sec. [5]). 

of 0, 5, 10, 20, a nd 30 km s " 1 . Further de tails can be found 
in lRvans et all (|2003l ) and lDufton et al.l (|2005l ). 

3.2. Limitations of the Be-type stellar analysis 

The nature of Be-type stars as discussed in the introduction 
are inconsistent with homogeneous, plane parallel, static 
model stellar atmospheres. The presence of a circumstellar- 
like disc, short-term variations and high rotation rates of 
possibly up to ~ 90% of the critical break-up velocity are 
characteristics that will affect the geometry and structure 
of the star but are not accounted for within the tlusty 
grid. For this reason the approach taken here is relatively 
simplistic, but a differential analysis should be sufficient to 
constrain the large range of nitrogen enrichments in our 
sample and provide some insights even in the light of these 
known issues. 

3.3. Stellar atmospheric parameters 

A static stellar atmosphere is characterised by the four 
parameters: effective temperature (T c ff), surface gravity 
(log g ) , microturbulent velocity (£ ) and chemical composi- 
tion. Due to the interdependence of these properties an it- 



1 .0/ 



1 .00 - 



0.98 - 



0.96 - 




0.94 1 1 1 

4530 4540 4550 4560 4570 4580 4590 



Fig. 1. Blending of Si in lines in NGC 346-008, which has a 
v sin i ~ 299 km s" 1 . Red lines indicate line centroids as cho- 
sen by the fitting procedure (for fixed relative wavelengths 
of the triplet). 
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erative procedure must be adopted when using the tlusty 
grids in order to obtain estimates for each parameter. For 
our targets, we used the appropriate SMC or LMC metallic- 
ity grid and initially assumed that the continuum contribu- 
tion from a circumstellar disc was negligible. We will sub- 
sequently refine our atmospheric parameters by attempting 
to allow for any such contribution. 

3.3.1. Effective temperature 

Effective temperatures were estimated using three different 
methods depending on the spectral quality, degree of ro- 
tational broadening and temperature regime. For objects 
with effective temperatures greater than 24 000 K, the He n 
absorption line at 4686A has a high temperature sensitivity 
and hence was adopted as our primary methodology. Other 
Hell lines at 4541 A and 4199A could also have been used, 
but in many of our objects neither of these lines were well 
observed. Therefore for consistency it was decided to only 
model the feature at 4686 A. The observed line profile was 
compared with theoretical spectra (corrected to allow for 
rotational broadening) from within the tlusty grids and 
the temperature estimated using a \ 2 technique. The qual- 
ity of the observed spectra typically leads to a typical un- 
certainty in T e ff of ±500 K but there were several objects for 
which larger uncertainties were appropriate due to the rel- 
atively poor quality of the agreement between observation 
and theory. As a note of caution, this method of determin- 
ing T c ff is dependent on the helium abundance of the object 
in question. The T e g values quoted for the objects in Tables 
[5] & El assume a normal helium abundance, but we can 
make an estimate to how the increase in helium abundance 
due to rotational mixing may affect the surface tempera- 
ture. Although He abundances have only been calibrated fo r 
slowly rotat i ng evol utionary models (jHeger fe Langerl2 000') . 
iBrott et all (|2011al) suggest a 10% enrichment of helium is 
possible for fast-rotating B-type stars. By using the tlusty 
grid of model atmospheres we find that an equivalent in- 
crease in the strength of the Hell 4686A absorption line 
leads to a temperature increase of <500K, hence we do not 
expect this to be a significant source of error. 

The silicon ionisation balance was used for cooler ob- 
jects without an observable He II spectrum, but where ab- 
sorption lines due to Sim (at 4552A, 4567A and 4574A) 
and cither Sin (at 4128A and 4130A) or Siiv (at 4116A) 
were observed. This methodology requires that the abun- 
dance estimates derived from the two ionisation stages 
agree. For objects whose spectra contained neither He II nor 
two silicon ionisation stages, effective temperatures were es- 
timated by requiring that the relative silicon to magnesium 
abundance was normal. The former was deduced from the 
Si ill multiplet at 4560A and the latter from the Mg II dou- 
blet at 448lA. The baseline abundances were taken from 
iHunter et all (|2007ft and are given in Table SJ Eleven ob- 
jects in the LMC used only this method, together with 
seven in the SMC. 

The use of three methodologies leads to the possibility 
of there being systematic discrepancies. Therefore where 
possible, we have compared for any given star the effective 
temperature estimates obtained using different methods. 
We note that this comparison was undertaken after the cor- 
rection discussed in Sect. l3.4l had been applied but our prin- 
cipal results would remain unchanged if we had used our 



Table 4. Present-day chemical com position of the 
Milky Way (iNieva fc Simon-Diad [2TJT1 . LMC and SMC 
(jHunter et al.l2007l) . Abundances are presented on the scale 
12 + \og[X/H]. 





MW 


LMC 


SMC 


c 


8.35 


7.75 


7.35 


N 


7.82 


6.90 


6.50 


O 


8.77 


8.35 


8.05 


Mg 


7.57 


7.05 


6.75 


Si 


7.50 


7.20 


6.80 



initial estimates. In the first instance, we adopted the other 
atmospheric parameters (log g and £) and estimate of disc 
contamination from our primary method and then redeter- 
mined only the effective temperature using the alternative 
methods. In Table [5J the effective temperature estimates 
are compared and in general, the agreement is reasonable 
particularly considering the caveats discussed in Sect. 13.21 
For six stars, the typical difference between the Hen and 
silicon ionisation balance was found to be 1 300 K with a 
standard deviation of 1 300 K. This is consistent with the 
He II, silicon ionisatio n equilibrium comparison performed 
bv lKilian etail (|l99lh . The comparison of the silicon ioni- 
sation balance and Si/Mg method was only possible in two 
cases but the values are in reasonable agreement. A compar- 
ison was inappropriate (and indeed normally not possible) 
between the Hen and Si/Mg methods, as the latter loses 
its temperature sensitivity above ~25 000 K. However there 
was one object, NGC 346-004, which could be analysed us- 
ing all three methodologies yielding estimates that had a 
range of 1 500 K. 

Following this comparison, a complete model atmo- 
sphere analysis (T c g , log g , £ and A [Disc] ) was re- 
performed for each star using the alternative methods for 
estimating the effective temperature. The results are again 
summarised in Table [SJ The difference in effective temper- 
ature estimates from the He n and Si equilibrium methods 
was now on average 1 700 K with a standard deviation of 
1100K. For the silicon ionisation equilibrium and Si/Mg 
methods the difference in one case was larger than in the 
previous comparison. For the one star, where all three meth- 
ods could be utilised, the range of effective temperature 
estimates was 1 600 K. 

For all three methods, observational uncertainties imply 
a typical error of ±500 K. However there are likely to be 
significant additional uncertainties due to the limitations 
in the theoretical methods (see Sect. 13.21) . which are very 
difficult to quantify. Where more than one method has been 
available, the estimates normally agree to within ±2 000 K 
implying that this might be a reasonable estimate of the 
effective temperature uncertainty. 

3.3.2. Surface gravity 

Surface gravities were initially estimated by comparing the 
observed Balmer lines, H 7 and H,5 with theoretical spectra 
from the tlusty grids; the estimates derived from both 
lines normally agreed to within 0.2 dex. Throughout our 
sample, the hydrogen line profiles tended to exhibit strong 
narrow absorption or em ission features at the line centre. 
iPorter fc Riviniusl (|2003l and references therein) discuss the 
structure of line profiles in Be-type stars in relation to the 
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Table 5. Comparison of atmospheric parameters and element abundances using different methods to estimate the 
effective temperature. 
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Notes. The method used is indicated as follows: Hell: profile fitting of the Hell absorption line at 4686A; Si: silicon ionization 
equilibrium, requiring at least 2 ionisation stages of Si; Si/Mg: uses the relative abundance estimates of silicon and magnesium 
see Sect. l3.3~Tl for more details. The estimated percentage contribution of disc emission (A[Disc]) to the continuum is also listed. 
The parameters from the primary methodology are listed first together with those from other effective temperature estimators, 
quoted as either 'Fixed' or 'Free'. 'Fixed' parameters are where other atmospheric parameters (log g and £) and estimates of the 
disc contamination are taken from the primary method, whereas 'Free' parameters allow all the parameters to change. ( -*- ) Assumed 
values of microturbulence 



kinematical nature of the circumstellar disc. The presence 
of a disc can result in single peaked or double peaked line 
profiles that may also be asymmetric and this was more 
prominent in the lower series members, H a - H 7 . Hence it 
was decided to adopt the estimates from the transition. 
The typical formal error in the fitting procedure was less 
than ±0.2 dex but there will be additional uncertainties due 
to the simplifications in our methodology. 

The use of different methods to derive the effective tem- 
perature (see Tablc[5]) led to a mean difference in the gravity 
estimates of 0.04 dex with a standard deviation of 0.14 dex 
for the six stars analysed by the Hen and silicon ionisa- 
tion equilibrium methods. Hence this would not appear to 
be a major source of uncertainty and we adopt a typical 
uncertainty in these estimates of ±0.2 dex. 

3.3.3. Microturbulence 

The microturbulent velocities were derived from the Si in 
multiplet at 4552 A, 4567Aand 4574A, by ensuring that the 
abundance estimates did not depend on line strength. For 
several objects, it was not possible to obtain a value of 
microturbulence that fulfilled this criterium and therefore 



an estimate of Okms -1 has been adopted, which minimises 
this depe ndency. This p henom enon has been discussed in 
detail by iHunter et al.l (|2007tl , where it was found in the 
normal B-type stellar population in the Magellanic Clouds. 
For some stars it was not possible to deduce a reliable mi- 
croturbulence as, for example, only one to two lines in the 
Si in multiplet could be measured, or the lines lay near to 
the linear part of the curve of growth. In these cases we 
ha ve adopted values th at are consistent with those found 
bv lHunter et al.l (|2009a|) for stars with similar gravities, viz. 
targets having a gravity between 3.5 - 4.5 dex were assigned 
^kms" 1 ; 3.0 - 3.5dex, £ = 10 kms" 1 ; 2.5 - 3.0dex, 
£ = 20 kms -1 . These cases are noted in Tables [S] and [5] 



The microturbulent velocities deduced using effec- 
tive temperatures derived from the He n profiles and 
the silicon ionisation equilibrium (see Table [5]) differ by 
4.0±5.0kms~ 1 , indicating the sensitivity of this quantity 
to the other atmospheric parameters. It is difficult to assess 
the reliability of these microturbulence estimates given the 
assumptions discussed in Sect. 13.21 However an uncertainty 
of the order of ±5 kms -1 would appear to be appropriate. 
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3.4. Corrections 

The initial estimates of the atmospheric parameters have 
been refined to allow for two effects, namely blending of the 
Mgn doublet with a nearby Aim doublet and continuum 
radiation from a circumstellar disc. 



3.4.1. Mgn line blending 

The procedure to estimate effective temperatures from the 
relative abundances of silicon and magnesium requires a re- 
liable measurement of the corresponding equivalent widths. 
The Mg II doublet is blended with an unresolved Al in dou- 
blet at 4479. 9A with this blending being more significant at 
higher effective temperatures. We have attempted to cor- 
rect for this blending as follows. FLAMES spectra of 30 
narrow lined B-type stars within our four clusters, cov- 
ering the effective temperature range ~ 12 000 - 33 000 K 
(jHunter et al.ll2007t iTThndle et al.ll2007[ ) were selected. For 
these objects, it was possible to independently measure the 
equivalent widths of Aim and Mgn features and then ob- 
tain a correction factor which was the ratio of the Aim 
equivalent width to the sum of the two equivalent widths. 

In Fig.[U this correction factor is plotted as a function of 
the effective temperature together with a polynomial least 
squares fit. This was then used to scale the observed Mgn 
equivalent widths (see Table [3]) of our sample and hence de- 
duce revised atmospheric parameter estimates for objects 
for which the Si/Mg method had been used. The effective 
temperature estimates were increased by typically ~ 200 K 
with a largest correction being ~ 700 K. Hence this correc- 
tion should not significantly increase the uncertainty in our 
eff ective temperature estim ates. 

iLvubimkov et al.l (|2005f ) performed a detailed analysis 
of the Alm/Mgii blending in a comparable sample of slow 
rotating Galactic B-type dwarfs. Comparison of o u r Fig . 
[2] and Fig. 2 from of the work of ILvubimkov et all ([2005) 
clearly illustrates that this blending effect is more promi- 
nent for the Milky Way sample than those in the LMC 
and SMC. As both samples are of comparable size (28 for 
Lyubimkov and 30 for this work) and span similar temper- 
ature ranges, it is reasonable to assume that the degree of 
blending is due to the metallicity environments. 



3.4.2. Correcting for disc contamination 

The issue of addressing the spectral contribution from a 
circumstellar disc around a Be-type star is not straight- 
forward. We have adopted a relatively simplistic approach 
that makes the following assumptions: 

— the contribution from the circumstellar disc is a feature- 
less continuum, that is a constant fraction of the stellar 
continuum over the spectral range considered. 

— being main sequence objects, our targets will have the 
same silicon abundance as normal B-type stars in the 
same cluster. 

Silicon abundances were estimated from the Sim triplet 
for all our targets and were compared with the baseline 
abundances (Table gj) for the SMC and LMC. The differ- 
ence (<5Si) was assumed to be due to the disc continuum 
and an iterative procedure was developed to allow for this. 
Initially, a fractional value for the disc continuum (C) was 
estimated from <5Si and the equivalent widths of all metal 
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Fig. 2. Effective temperature is plotted against a cor- 
rection factor, EW(Allll)/EW(Alm+Mgll), derived from 
the Aim and Mgn absorption lines at 4479. 9A and 
448 1.0 A respectively. This correction factor was obtained 
from the analysis fo r a s ample of 30 B-type sta rs taken from 
iHunter et all (|2007h and lTrundle et al.l (|2007h and based on 
the polynomial fit given by the solid black line. 



lines were then scaled by the reciprocal of (1 — C). These 
revised estimates of the stellar equivalent widths were used 
to estimate the effective temperature for stars using the sil- 
icon ionisation balance or Si/Mg methods. For stars with 
temperatures estimated using the He 11 spectrum, the line 
profiles were corrected by subtracting the estimated frac- 
tional contribution and renormalising the feature. A simi- 
lar procedure was applied to the H,5 lines and the surface 
gravity re-estimated. To complete the iteration, microtur- 
bulences were re-calculated using the corrected silicon line 
strengths. With these new atmospheric parameters, silicon 
abundances were re-estimated, resulting in a revised disc 
contribution. This iteration procedure was repeated until 
the values <5Si were within ±0.1 dex. After applying these 
corrections to those targets where the disc contribution was 
non-zero, the typical increase in effective temperature was 
~ 1 600AT, the increase in surface gravity was ~ 0.45 dex 
and the increase in microturbulence was < 5kms _1 . 



3.4.3. Photospheric abundances 

Table [5] summarises the photospheric abundance estimates 
(or upper limits) for our targets with all values rounded 
to the nearest 0.1 dex. Our principal aim has been to try 
and constrain the evolutionary status of our Be-type sam- 
ple and hence their nitrogen abundances are crucial. The 
methodology has also led to estimates of the silicon and 
in some cases magnesium abundances. However these are 
not expected to be significantly affected by mixing from a 
hydrogen burning stellar core. Additionally the spectra of 
these elements were used to deduce atmospheric parameters 
and the amount of continuum emission from a circumstel- 
lar disc and hence these abundance estimates cannot be 
cons idered to be indep endent. 

iFraser et aD (|2010h used the Nil line at 3995A for a 
sample of supergiant stars and investigated how their abun- 
dance estimates differed from those deduced using all avail- 
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able Nil lines. They concluded that a conservative esti- 
mate of the random error was ±0.15 dex. Additionally for 
six stars with effective temperatures between 13 000 and 
25 000 K, they found that the difference in the abundance 
estimate from the line at 3995 A and the other N n lines was 
0.03±0.08 dex. Hence we do not expect that the use of only 
one line should lead to significant systematic errors. As an 
additional test, the Be-type sample was searched for addi- 
tional Nil lines, with identifications being made in 11 stars. 
The other lines observed were 4630A, although this was of- 
ten contaminated by an emission feature around 4629A, 
and 4447A; in two cases a line at 4601 A was also visible. 
Including the abundance estimates from these features for 
objects with T c g <26000 K led to nitrogen abundances that 
were in general consistent with those estimated from the 
3995Aline alone, apart from NGC 346-004, where it was 
found to be 0.2 dex higher. However as discussed in Sect. 
14.21 this star exhibits shell like features thereby limiting the 
reliability of these measurements. 

There could also be systematic errors due to, for exam- 
ple, errors in the atmosp heric paramet e rs. Th ese have been 
extensively discussed by iHunter et al.l (|2007) and vary de- 
pending on the strength of the feature and on the adopted 
atmospheric parameters. Typically such errors are of the or- 
der of 0.1 - 0.2 dex (but on occasions can be larger), which 
is consistent with the range of abundance estimates listed 
in Table [5] for individual stars. 

Combined these two sources of error would typically 
lead to an uncertainty of ±0.2 dex. However there may be 
additional systematic errors due to for example the simpli- 
fications made in estimating the continuum contamination 
due to the disc. These are difficult to quantify but must 
be borne in mind when interpreting these results. For ex- 
ample, several of our targets appear to have negligible disc 
contamination despite their Be nature. In Fig. [3] we show 
the effect that this disk contamination has on our derived 
nitrogen abundances and one can see an approximate cor- 
relation such that every error of 10% in the disc contami- 
nation results in an error of 0.1 dex in nitrogen abundance. 
Therefore for some objects our uncertainties may be as high 
as 0.4 dex. 

4. Discussion of results 

4.1. Atmospheric parameters 

Atmospheric parameters and photospheric abundances are 
presented for a total of 30 stars, 16 from the SMC and 
14 from the LMC in Table ED The distribution of effective 
temperatures and gravities are shown in Fig. U and indi- 
cate that the SMC and LMC targets span similar effective 
temperature ranges. The surface gravities and microturbu- 
lent velocities also span similar ranges, although for the 
latter this may be due to the significant proportion (43%) 
of adopted values (see Sect. 13.3.3)) . 

4.2. Nitrogen abundances 

For the SMC sample, nitrogen abundance estimates have 
been obtained for Be stars with relatively high projected 
rotational velocities (greater than 2 00 km s -1 ), in contras t 
to the results of the B star sample of lHunter et al.l (|2009a|K 
which normally had only upper limits. Inspection of Table 
EO also reveals that these Be stars tend to have surface grav- 
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Fig. 3. The change in nitrogen abundance after applying 
the contamination factor is plotted against the estimated 
continuum contribution from the presence of a circumstel- 
lar disc . Both LMC and SMC stars are included (circles). 
Stars showing shell-like features are noted as open circles. 
A linear least squares fit (solid line) has been performed on 
all points and indicates a good correlation. 
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Fig. 4. Effective temperature plotted against surface grav- 
ity estimates for all stars listed in Table E3 Points are plot- 
ted for combined cluster samples; circles for the SMC and 
diamonds for the LMC. Open symbols represent stars with 
shell-like features. 



ities lower than one would expect for (near) main sequence 
stars. 

Therefore, we have reviewed the spectra of our Be- 
type sample and find the presence of previously over-looked 
shell-like Balmer lines in eight stars; NGC 346-004, 008, 024, 
036 and NGC 2004-027, 035, 067, 096. Since we are presum- 
ably observing these stars through their equatorial discs it 
is likely that our simple assumption that we observing the 
sum of the photospheric spectrum and disc continuum is no 
longer valid, leading to an underestimate of the stellar grav- 
ities. Although it would be possible to correct these surface 
gravities, the impact of absorption by the disc on the ni- 
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Table 6. Estimates of the atmospheric parameters, disk contamination and silicon, magnesium, nitrogen abundances for 
the Be-type sample. 
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trogen spectrum is largely unknown. Therefore we have not 
sought to correct the nitrogen abundances for these shell- 
stars but flag them in Fig.s[3]&|3]as they may be unreliable. 
A further eight stars may also show signs of disc absorption 
(NGC 346-009, 023, NGC 330-025, 031, NGC 2004-023, 056, 
N 11-073, 074) and their results should be regarded with 
some caution. 

An apparent correlation between effective temperature 
and surface gravity is seen in Fig. |H partially due to 8 ob- 
jects with log g=2.5- 3.4 and with T c g <20 000 K, namely 
NGC 346-008, NGC 346-024, NGC 330-025, NGC 330-069, 
NGC 2004-025, NGC 2004-048, NGC 2004-039, NGC 2004- 
089. The stellar parameters of these objects are consistent 
with objects towards the end (or beyond) their main se- 
quence lifetime and might therefore show a stronger enrich- 
ment of Nitrogen. As discussed in Sect. l3.4.Ul there is a pos- 
sibility that use of only the single N II 3995A line could re- 
sult in an underestimation of the true nitrogen abundance. 
Therefore these objects were again checked for additional 
N II lines. Out of the 8 objects, 4 showed at least one extra 
identifiable N II line, and by taking the averaged abundance 
from all measured N II lines 3 objects showed an enhance- 
ment to their nitrogen abundance of approximately 0.1 dex. 



The remaining object is noted as having shell like features 
and hence the estimates may be unreliable. 

The nitrogen abundance estimates in both metallicity 
regimes span significant ranges: SMC: 6.7 - 7.8 dex, LMC: 
7.0 - 8.0 dex. The largest enrichments in both metallic- 
ity regimes are of a similar magnitude, being ~1.3 and 
-1.1 dex, for the SMC and LMC respectively. For the SMC 
eight objects showed no evidence of the Nil line at 3995A 
and the upper limits to their absolute nitrogen enhance- 
ments (0.2 - 1.2 dex) indicate that at least in some cases the 
nitrogen enrich ment must be small. This is consistent with 
an analysis by iLennon et al.l ([2005) who found effectively 
no nitrogen enrichment for a pole-on Be-type star in the 
SMC cluster NGC 330. Only one object (NGC 2004-092) 
in the LMC yielded a nitrogen upper limit, which implied 
effectively no enrichment (<0.1dex). 

A surprising feature of our results is that the nitrogen 
abundances of normal B and our Be-type stars are found 
to have similar values. The normal B-type stars are defined 
as those which exclude Group 1, N-normal fast rotators, 
and Group 2 , N-rich slow rotating stars, (as discussed by 
iHunter et all l2009al. . In Fig. El the range of B-type and 
Be-type nitrogen abundances are given for the SMC and 
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LMC. For the SMC, there are a significant number of up- 
per limits applied to the nitrogen measurements for fast 
rotating stars in both samples. Nevertheless, this figure il- 
lustrates that even for the SMC, Be-type and B-type stars 
occupy a similar range of nitrogen abundances, and hint 
at a stronger enhancement at lower metallicity. This com- 
parison between Be- and B-type stellar abundances raises 
an important question; If Be-type stars are all fast rotators 
(70% of critical velocity. IPorter fc Riviniusl 120031 and ref- 
erences therein) then why are they not more N-rich than is 
observed? We return to this in Sect. [SJ 

As discussed previously one would consider more 
evolved giant and supergiant objects to exhibit a 
stronger enhanc e ment to their nitrogen abundance. 
iLvubimkov et all (|2011l) observed such nitrogen enhance- 
ments in A- and F-type supergiants which could have ob- 
tained their nitrogen enrichment from rotational mixing 
during a B-type MS lifetime. To examine this further it 
would have been useful to investigate the relationship be- 
tween observed nitrogen enrichment and fundamental pa- 
rameters such as the mass and main sequence lifetime of 
our objects. For Be-type stars, however, obtaining the in- 
trinsic luminosity is not a straight forward task, due to the 
presence of cir cumstellar material and effects such as grav- 
ity darkening (jTownsend et al.l [20041) . and as a result we 
were unable to obtain any useful constraints. 

If we consider the relative nitrogen enrichments in ref- 
erence to the cluster ages, NGC346 & the Nil region 
are of roughly the same age of a few million years, whilst 
NGC 2004 & NGC 330 are older but of comparable ages. 
However we find there to be no correlation between the 
abundances e stima tes and cluster ages. As discussed by 
lEvans et al.l (|2006lb the majority of the targets in the 
FLAMES survey are probably field stars, whose ages may 
well be different from that of the cluster. 

4.3. Projected rotational velocities 

Hun ter et al.l (|2008l ) have previously discussed the pro- 
jected rotational velocities for their FLAMES sample, 
including the Be sample considered here. In particular 
they compared their result s with a previous study of 
IMartavan etail (|2006l I2007D of a large sample of B- and 
Be-type stars in NGC 2004, and NGC 330. 

We have modelled our pr ojected rotational velocities 
and those of Martav an et al.l to infer the intrinsic rota- 
tional velocity distrib utions, which w ere assum ed to follow 
a Gaussian profile - see lHunter et al.l fsee l2008h for details of 
the m ethodology. Fig. [5] illustrates the results for the sam- 
ple of IMartavan et all the quality of the fit for our sample 
was somewhat worse, due to our smaller sample sizes. The 
inferred rotational velocity distributions are summarised in 
Table [7] and imply that the majority of these Be-type stars 
are currently rotating at velocities that are considerably 
lower than their cr i tical ve locities. From evolutionary mod- 
els of iBrott et al.l (|2011al ) the critical velocity is typically 
600 - 650 km s -1 for our mass range, which would suggest 
our objects have velocities <50% of critical. We return to 
these results in Sect. [5] 

5. Comparison with evolutionary models 

While the evolutionary status of Be-type stars is still un- 
certain it is interesting to compare our derived nitrogen 



Table 7. Mean rotational velocities, - widths and stan- 
dard deviations (cr) for an assumed Gaussian distribution 
of rotational velocities. 



Sample n wo - cr 

(kms" 1 ) (kmV 1 ) (kms" 1 ) 

This work: LMC 19 210 100 70 

This work: SMC 46 260 140 100 

Mar07: NGC 2004 47 190 205 145 

Mar06:NGC330 131 310 190 135 



Notes. Velocitie s are considered for b oth the LMC and SMC 
by this work fsee | Hunter et~ al. 2009b, for further details) and 
that of lMartavan et all Giventhe uncertainties associated with 
these values they have all b een rounded to the nearest 5 km s" 1 . 
Mar06 and Mar07 reference IMartavan et all ||2006|.|2007T ) respec- 
tively. 

abundances with the predictions of stellar evolution mod- 
els in the context of rotational mixing. If one assumes that 
all Be-type stars are fast rotators, but are otherwise nor- 
mal B-type stars, then clearly they should be important 
tests of the efficiency of rotational mixing. Here we com- 
pare our abundance estimates with those predicted from a 
simulated population of B-type stars base d on theoretical 
evolutionary models of iBrott et al.l (|201 lai ibh . 

5.1. Model assumptions 

The code, STARMAKER, (|Brott et al.l 1201 lbf ) simulates 
populations of early-type sing le stars, based on t he evo- 
lutionary models presented in IBrott et al.l (|2011aj) , which 
take into account effects of rotation and mass loss. 
STARMAKER predicts physical quantities (e.g. atmo- 
spheric parameters and abundances) for a prescribed sam- 
ple of typically 1-2 million of stars with initial mass, 
initial rotational velocity and age drawn from the corre- 
sponding distribution function (Salpeter IMF, initial distri- 
bution of rotational velocities, star formation history (SFH) 
and random inclination angl e). We adop t ed the parame- 
ters used in the simulations of IBrott et al.l (|2011br i as these 
should be appropriate to the current VLT-FLAMES sam- 
ple. In particular the underlying distribution of initial stel- 
lar rotational velocities reproduced the observed distribu- 
tion of projected rotational velociti es of the VLT-FLAMES 
B-sta rs. As discussed in Sect. 14.31 IMartavan et al.l ([2006, 
120071 ) also observed a sample of B- and Be-type stars in 
NGC 2004, and NGC 330. We note that the use of this dis- 
tribution with our simulations do not change our results. 

The use of observed current projected rotational veloc- 
ities to simulate the initial distribution of rotational veloc- 
ities assumes that there is a relatively small change during 
the main s equence lifetimes. Th is is confirmed by the sim- 
ulations of IBrott et"al~l (|2011bD and illustrated in Fig. [71 
which shows effectively no change in its rotational velocity 
during the main sequence lifetime for a twenty solar mass 
star. 

5.2. Predictions 

Selection effects were applied to the simulated popula- 
tion in order to represent the observed sample. Typically 
^75% of the simulated sample were removed to match the 
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observed faint magnitude cut-off (i.e. exclusion of lower 
mass- stars), in addition to stars with ages beyond their 
main sequence lifetime and those with T e ff>35 000K and 
log g <3.2 d ex. Other po s sible selecti on effects h a ve bee n 
discussed in lEvans et al.l (|2006t ) and iBrott et all (|2011bft . 
However given the relatively simple nature of our compar- 
isons we do not believe that these will be important. 

With the assumption the population synthesis mod- 
els should be able to reproduce our observed objects, we 
first investigate whether Be-type stars are a population of 
rapidly rotating B-type stars. In Fig.s|8]&[9]simulations are 
presented assuming that Be-type stars rotate around 70% 
or 40% of their critical rotational velocity. The simulations 
with 70% clearly predict too high nitrogen abundances for 
the LMC, while for the SMC the picture is more mixed with 
some stars being well matched to the range of predicted 
abundances, while others are too low. When we consider 
the simulations for stars at around 40% of their critical 
velocity the LMC gives a better agreement, but again the 
SMC picture remains mixed and as already noted is com- 
plicated by the presence of more upper limits in the SMC. 

Based on the nitrogen abundances, and assuming our 
evolutionary models are realistic, all Be-type stars do not 
rotate close to critical velocity as is often assumed. This 
is consistent with the mean velocities found for different 
Be-type samp l es in Table [7] Note that as discussed by 
iMokiem et al.l (|2006| ). observed projected rotational veloc- 
ities do not well constrain details of the distribution of the 
underlying rotational velocities. However they do constrain 
the mean rotational velocity and hence the values in Table 
[7] should be reliable. 

Although there are significant uncertainties in our esti- 
mated nitrogen abundances as discussed in Sect. 13.4.31 we 
believe that they are not compatible with stars that have 
spent most of their evolutionary history rotating at close to 
their critical velocity. As discussed above, simulations for 
such cases imply that effectively all the sample should have 
large nitrogen enhancements, which is not observed. We 
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Fig. 8. Estimated nitrogen abundances and projected ro- 
tational velocities plotted for the SMC (top panel) and 
LMC (bottom panel) Be-type stellar samples before and af- 
ter contamination correction (blue and black respectively) . 
Simulations are based on the rotational velocity distribu- 
tion of Hunter et al. with only B-type stars rotating with 
critical velocity between 65% and 75% as to represent a 
sample of stars rotating around 70% of critical velocity. 
Marked on each plot as a black dashed line is the baseline 
nitrogen abundance. Simulations predict too high nitrogen 
abundances for the LMC, with a mixed view given for the 
SMC. 
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can postulate different scenarios to explain our estimated 
nitrogen abundances, including: 

— On average Be-type stars rotate faster than normal B- 
type stars but have velocities significantly less than crit- 
ical. 

— If Be-type stars have had rotational velocities close to 
critical values, this phase must have been for a relatively 
short period of their main sequence lifetime. 

For the first hypothesis, we are faced with the fun- 
damental difficulty of why there are two types of stars. 
Certainly on average Be-type stars may have higher ro- 
tational velocities than B-type stars. However the range 
of velocities for Be-type stars implied by the a- va lues in 
Table [7] and for B-type stars in lHunter et al.l (|2009af ) imply 
that there will be a large overlap in the range of rotational 
velocities of the two groups. In turn this will imply that a 
substantial fraction of Be-type stars have current rotational 
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Fig. 9. Estimated nitrogen abundances and projected ro- 
tational velocities plotted for the SMC (top panel) and 
LMC (bottom panel) Be-type stellar samples before and af- 
ter contamination correction (blue and black respectively). 
Simulations are based on the rotational velocity distribu- 
tion of Hunter et al. with only B-type stars rotating with 
critical velocity between 35% and 45% as to represent a 
sample of stars rotating around 40% of critical velocity. 
Marked on each plot as a black dashed line is the base- 
line nitrogen abundance. Simulations predict a better agree- 
ment with LMC abundances and again a mixed view given 
for the SMC. 



velocities that are less than normal B-type star mean ve- 
locity. 

If we assume that Be-type stars have experienced near 
critical rotation, we can use our grid of evolutionary mod- 
els to infer an upper limit for the fraction of the main se- 
quence lifetime. In Fig. [JJ the evolution with time is shown 
for a 20M Q SMC B-type star with an initial rotational ve- 
locity ~ 400 km s -1 . The nitrogen enrichment occurs rela- 
tively early in the main sequence lifetime, the majority in 
the first two million years (corresponding to less than one 
quarter of the main sequence lifetime). At this stage the 
nitrogen enhancement is more than 1.0 dex corresponding 
to the largest enhancements found in our sample. Therefore 
if we take the average nitrogen enhancements found in our 
sample and we believe the theory of rotational mixing to be 
correct the star could only exist close to critical velocity for 
within the first ten percent of the main sequence lifetime. 
Consideration of oth er models within our gri d or those of 
other investigations (jMevnet fc Maederl I2002T ) imply simi- 



lar results with the nitrogen enrichment (and hence a short 
period of close to critical velocity) occurring early in the 
main sequence phase. Alternatively, it is possible that the 
stars have undergone a recent spin up event, possibly due 
to an interacting binary system and that large-scale nitro- 
gen enrichment (by rotational mixing) has not yet taken 
place. This would require inefficient mass transfer as oth- 
erwise this could contribute to the gaining star's nitrogen 
abundance. 

Recent work by iHuang et al.l (|2010h discusses the rota- 
tion of Galactic B-type stars and infer a threshold for Be- 
type star phenomena by considering the maximum value 
of the rotational velocity to the critical velocity for their 
normal B-type star population. This is found to vary with 
stellar mass from ~ 95% for stellar masses < 4M Q to ~ 
60% for the higher masses (> 8.6 M Q ) appropriate to our 
sample. Although we cannot explicitly determine at what 
rate of critical velocity our sample appears to be rotating 
at, from our population synthesis mode ls and measur ed pa- 
rameters the fraction of v cr i t found by IHuang et al.l would 
be consistent with Be stars in our mass range." 



i 6. Conclusions 

% We have presented projected rotational velocities, atmo- 
| spheric parameters and photospheric abundance estimates 
z of three elements, including nitrogen for a small sample 
of Be-type stars in the SMC and LMC in order to com- 
pare with current evolutionary models. Our principal con- 
clusions are: 



— The estimated nitrogen enrichments are found to be 
similar for both Be- and B-type stars. 

— Nitrogen abundances and projected rotational velocities 
have been c ompared with popu lation synthesis models 
produced bv lBrott et al.l (|2011bf ). They are inconsistent 
with a population that have spent most of their lifetime 
at near critical rotation velocities. 

— Be-type stars could have spent their main sequence life- 
time rotating at velocities on average larger than B 
stars but still significantly less than the critical velocity. 
However the spread in velocities observed then implies 
that a substantial proportion of Be-type stars are cur- 
rently rotating at velocities smaller than B stars. 

— Given the nitrogen enrichments found in our samples, 
Be stars could only experience near critical rotational 
velocities for a relatively short period of their main se- 
quence lifetime (less than 10%). 
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